Strongly paired fermions: Cold atoms and neutron matter 
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Experiments with cold Fermi atoms can be tuned to probe strongly interacting fluids that are 
very similar to the low-density neutron matter found in the crusts of neutron stars. In contrast to 
traditional superfluids and superconductors, matter in this regime is very strongly paired, with gaps 
of the order of the Fermi energy. We compute the T — equation of state and pairing gap for cold 
atoms and low-density neutron matter as a function of the Fermi momentum times the scattering 
length. Results of quantum Monte Carlo calculations show that the equations of state are very 
similar. The neutron matter pairing gap at low densities is found to be very large but, except at 
the smallest densities, significantly suppressed relative to cold atoms because of the finite eflFective 
range in the neutron-neutron interaction. 
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Strongly paired fermions are important in many con- 
texts: cold Fermi atom experiments, low-density neutron 
matter, and QCD at the very high baryon densities po- 
tentially found in the center of massive neutron stars. De- 
veloping a quantitative understanding of strongly paired 
Fermi systems is important since they offer a unique 
regime for quantum many-body physics, relevant in very 
different physical settings including the structure and 
cooling of neutron stars. Constraining neutron matter 
properties can also be important in understanding the 
exterior of neutron-rich nuclei by constraining parame- 
ters of nuclear density functionals. 

Cold-atom experiments can provide direct tests of the 
equation of state and the pairing gap in the strongly 
paired regime, and hence provide a crucial benchmark 
of many-body theories in these systems. We consider a 
system of two fermion species and a simple Hamiltonian 
of the form 



H 



2m 



(1) 



where i and j represent spin up and down particles, re- 
spectively. In cold atoms the interaction v{r) can be 
tuned through Feshbach resonances to be very attrac- 
tive, and to produce a specific scattering length. By 
varying the scattering length one can sweep through 
various values of Fermi momentum kp times scatter- 
ing length a from the BCS side (small —kpo) to uni- 
tarity {—kpa = oo) and beyond. In many experiments 
^Li atoms are used, and the effective range r^. between 
the atoms is nearly zero. A variety of beautiful experi- 
ments have been performed recently on cold atoms in the 
strongly paired regime [DiaialllllSllSllIllIllIllIO]. 

The neutron-neutron interaction, in contrast, is gener- 
ally quite complicated, with one-pion exchange at large 
distances, intermediate range spin-dependent attraction 
dominated by two-pion exchange, and a short-range re- 
pulsion. At very low densities, though, as found in neu- 



tron star crusts or the exterior of neutron-rich nuclei, the 
scattering length and effective range are most crucial to 
the physical properties of the system. The presence of a 
short-range repulsive core is important primarily in that 
it prevents collapse to a higher-density state. 

A strongly attractive zero-range interaction was pro- 
posed as a simple model of neutron matter even before 
the recent remarkable cold atom experiments |llj . In 
low-density neutron matter the scattering length is very 
large, « —18.5 fm, much larger than the typical sep- 
aration between neutron pairs. The effective range is 
much smaller than the scattering length, re « 2.7 fm, so 
\re/a\ ~ 0.15, but only at very low densities is the effec- 
tive range much smaller than the interparticle spacing. 

To the extent that the effects of finite range in the in- 
teraction can be neglected, cold atoms and neutron mat- 
ter are 'universal' in the sense that the properties of the 
system depend only upon the product of the Fermi mo- 
mentum and the scattering length. Experiments have 
been performed that probe the sound velocity |T] and 
collective excitations superfluidity [3 and critical 
temperature |4j, phase separation and phase diagram 
[3 [am HE] and RF response [10]. 

We have performed fixed-node quantum Monte Carlo 
(QMC) calculations for both cold atoms and neutron 
matter. In each case, the trial wave function is taken 
to be of the Jastrow-BCS form with fixed particle num- 
ber and periodic boundary conditions: 



i<3 



(2) 



The BCS pairing function (j){r) is parametrized with a 
short- and long-range part as in 12j. Since the ground- 
state energy in a fixed-node calculation is an upper bound 
to the true ground state energy, we optimize the parame- 
ters to obtain the lowest fixed-node ground-state energy, 
as in [T^ . 
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FIG. 1: Zero-temperature equation of state for cold atoms 
and neutron matter. Near zero density we show the analytical 
expansion of the ground-state energy of a normal fluid, and 
at high density we show the cold atom result at unitarity 
{kpa = oo, arrow). QMC calculations are shown as circles 
and squares for neutron matter and cold atoms, respectively. 



The interaction for cold atoms is taken as v{r) 



-Vq 



2h^ M 

m cosh^(/^r) 



with jj, = 24:/ tq, or an effective range 

of ro/12, with 1/p = (4/3)7rrQ. The interaction range 
is small enough not to significantly affect the energy 
or pairing gap from the BCS regime to unitarity. For 
the neutron-neutron interaction, we take the s-wave 
part of the AV18 [13 interaction. This interaction fits 
nucleon-nucleon scattering very well at both low- and 
high-energies. For our purposes the important thing is 
that the scattering length and effective range are cor- 
rectly described. We use this interaction only between 
spin up-down pairs, which sets the interaction in the 
L = = ±1 pairs to zero. We correct for the artifi- 
cial attraction in the L = 1 , M = pairs perturbatively. 
This correction is 10% for the ground-state energy at the 
largest densities considered, and typically much smaller. 
The correction to the pairing gap is always smaller than 
the statistical error in the calculation. 

The T = equations of state for cold atoms and neu- 
tron matter are compared in Fig. [l] The horizontal axis 
is kpa, with the equivalent Fermi momentum hp for neu- 
tron matter shown along the top. The vertical axis is the 
ratio of the ground-state energy to the free Fermi gas 
energy (Epc) at the same density; it must go to one at 
very low densities and decrease as the density increases 
and the interactions become important. The curve at 
lower densities shows the analytical result |14j for normal 
matter: E/Epc = ^akp+ ^ (11 -2 In 2) {akpf. 
This curve should be valid at very low densities. While 
it ignores the contributions of superfluidity, these are ex- 



ponentially small in (\/kpa). 

The neutron matter and cold atom equations of state 
are very similar even for densities where the effective 
range is comparable to the interparticle spacing. Hence 
cold atom experiments can tell us something rather di- 
rectly about the neutron matter equation of state. Near 
kpa = —10 the energy per particle is not too far from 
QMC calculations (THIH] and measurements 16J of the 
ratio ^ of the unitary gas energy to Epq; previous calcu- 
lations give ^ = 0.42(1). Extrapolations of recent QMC 
calculations to re = and also AFMC calculations sug- 
gest that C = 0.40(1) [17] (arrow in Fig. [l}. 

The results near kpa = — 10 for neutron matter are 
compatible with previous calculations of the neutron 
matter equation of state at somewhat higher densities 
{kp > 1 fm-i) [TSllTnillD]. Results shown are for 66 par- 
ticles in periodic boundary conditions; calculations have 
also been performed near iV = 20, 44, and 90. Based on 
these results, finite-size effects for TV = 66 and beyond 
are expected to be quite small, of the order of a couple 
percent. Calculations of the cold atom equation of state 
are very similar to those reported previously in [15j and 
PT| : the energies reported here are slightly lower (up to 
« 10% in some cases) because of larger system sizes and 
better optimizations. 

Realistic microscopic calculations that incorporate 
strong pairing thus provide important constraints on the 
neutron matter equation of state in the subnuclear satu- 
ration density regime. Skyrme models or more generally 
density functionals are used, for example, to determine 
the structure of neutron star crusts i22j and the neu- 
tron skin thickness of nuclei A realistic treatment 
of these problems should incorporate the physics of the 
rapid transition of neutron matter from nearly free par- 
ticles to a strongly paired system at very low densities. 

The pairing gap is the other fundamental zero- 
temperature property of superfluid systems. Calcula- 
tions of the s-wave pairing gap in neutron matter have 
varied enormously over the past 20 years [531 123 • The 
difficulties in accurately calculating corrections to the 
BCS pairing gaps in the strongly paired regime are 
significant, and hence calculations of the pairing gap 
[H 1151 HSl 123 IMIIMl [3D] can differ by large factors 
(from 4 to 10) in the low-density regime. Cold atom ex- 
periments can provide a critical test of theories of the 
pairing gap in this regime. We first compare our cal- 
culations of the pairing gap in cold atoms and neutron 
matter, and then compare with previous results. 

We calculate the pairing gap from the odd-even energy 
staggering A = E{N+1)-{eIn)+E{N+2)) /2, where N 
is an even number of particles. Finite-size effects for the 
pairing gap are considerably larger than for the ground 
state energy. In order to estimate the convergence of the 
gap to the continuum value with increasing N we have 




FIG. 2: Neutron matter pairing gap at kpa — —10 versus 
particle number in periodic boundary conditions, BCS and 
QMC calculations. 



FIG. 3: Superfluid pairing gap versus kpa for cold atoms 
(r-e ~ 0) and neutron matter (|re/a| « 0.15). BCS (solid 
lines) and QMC results (points) are shown. 



solved the BCS equations: 



A(k) = -V(k|F|k')— ^^^'^ 



e(k) 



Ve(k)2 + A(k)2 



(3) 



in periodic boundary conditions for different (N). 

The hne in Fig. |2] is the continuum BCS result for 
kpa = —10, and the open symbols are the solutions of the 
BCS equations for different {N). The continuum results 
are nearly identical for the AV18 interaction and the sim- 
ple cosh potential adjusted to yield the same scattering 
length and effective range. For the finite systems BCS re- 
sults are shown for the cosh potential. Unlike the case of 
cold atoms near unitarity, where —kpa » 1 and ~ 0, 
the BCS gap shows sizable oscillations for small numbers 
of particles. The BCS value approaches the continuum 
limit (straight line) near = 66, and oscillations from 
that point on are fairly small, comparable in size to the 
statistical error in the QMC calculations. We also show 
as solid points the gaps obtained from particle-projected 
BCS wave functions in variational Monte Carlo calcula- 
tions and the odd-even staggering formula. The projec- 
tion to definite particle number is a small effect. 

The lower points in Fig. [2] are QMC results for 
kpa — —10. At very small values of N the gap is quite 
large, as is also seen in the BCS calculations. This is 
due to the coarse description of the Fermi surface in such 
small systems; the momentum grid spacing in occupied 
states is similar in magnitude to the Fermi momentum. 
When the pairing is very strong, as in cold atoms in the 
unitary regime, this coarse description is not too critical. 



However for weaker coupling or the larger effective range 
in neutron matter this becomes more important. The 
gap in both BCS and QMC calculations reaches a mini- 
mum near 44 particles (near the midpoint between closed 
shells at 38 and 54) , and then increases to values near the 
continuum limit. Pairing gap results for A^ = 66 — 92 are 
consistent within the statistical errors. 

For all values of A^ the gap is considerably smaller than 
the BCS results. For comparison, at unitarity in cold 
atoms BCS calculations give a gap of 0.69 Ep while the 
QMC result is 0.50(3) Ep.'ST] These calculations are in 
good agreement with recent polarized cold atom exper- 
iments [HI [32] • For cold atoms the BCS equations will 
produce the exact gap in the BEC limit where the pairs 
are strongly bound. No such limit is relevant for a finite- 
range interaction. 

In Fig. |3]we plot the pairing gap as a function of kpa 
for both cold atoms and neutron matter. BCS calcula- 
tions are shown as solid lines, and QMC results are shown 
as points with error bars. QMC pairing gaps are shown 
from calculations of A^ = 66 — 68 particles. For cold 
atoms away from unitarity the pairing gaps are smaller 
than calculated previously [2r, due to more complete op- 
timizations and because these larger simulations reduce 
the finite-size effects. 

For very weak coupling, —kpa « 1, the pairing 
gap is expected to be reduced from the BCS value by 
the polarization corrections calculated by Gorkov [SS] 
A/Abcs — (l/4e)^/'^. Because of finite-size effects, it 
is difficult to calculate pairing gaps using QMC in the 
weak coupling regime. The QMC calculations at the low- 
est density, kpa = —1, are roughly consistent with this 
reduction from the BCS value. At slightly larger yet still 
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very strongly paired, and both have pairing gaps of the 
order of the Fermi energy. Experiments on the cold-atom 
equation of state would be very valuable in constraining 
the neutron matter equation of state. Pairing gaps in 
neutron matter are found to be suppressed compared to 
cold atoms and BCS theory, but much larger than in 
most other approaches. Again, cold-atom experiments 
could provide very valuable tests of many-body theories. 
It could be very important to explore finite-range effects 
experimentally using other atomic systems. 

Note added in proof. We recently became aware of new 
calculations of the equation of state and pairing gap for 
cold atoms using auxiliary field Quantum Monte Carlo 
techniques [35) . Their results are similar to, but slightly 
different than, those presented here. 



FIG. 4: Superfluid pairing gap versus kpa for neutron matter 
compared to previous results. 



small densities, where —kpa — 0(1) but fc^^e << 1 for 
neutron matter, one would expect the pairing gap to be 
similar for cold atoms and neutron matter. The results 
at kpa = —2.5, where kpr^ ~ 0.35, support this expec- 
tation. Beyond that density the effective range becomes 
important and the QMC results are significantly reduced 
in relation to the cold atoms where ~ 0- 

These results for the pairing gap are compared to 
selected previous results in Fig. |4] The results of 
our calculations are much larger than the diagrammatic 
[?7l EH] and renormalization group [29 approaches. 
As these approaches assume a well-defined Fermi surface 
or calculate polarization corrections based upon single- 
particle excitations it is not clear how well they can de- 
scribe neutron matter in the strongly paired regime, or 
the similar pairing found in cold atoms. 

The results here are significantly smaller than the 
AFDMC results of Fabrocini et al. [3D]. These calcu- 
lations are somewhat similar to those reported here. The 
disadvantage of the AFDMC approach is that it does 
not provide a variational bound to the energy, and hence 
the wave functions are chosen from another approach. 
In the calculation of Ref. 30J the wave function was 
taken from a correlated basis function approach that in- 
cluded a BCS initial state. The pairing in that variational 
state is unusually large, and in fact increases as a frac- 
tion of Ep when the density is lowered. The advantage 
of the AFDMC approach is that it includes the full in- 
teraction rather than the simple s-wave interaction used 
here. AFDMC calculations with larger particle numbers 
are underway. |34j 

In summary, we have calculated the T = equations 
of state and pairing gaps for cold atoms and neutron 
matter. These systems are quite similar in that both are 
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